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Abstract 

The  transonic  canard-wing  Interaction  problem 
is  simulated  using  modified  small  disturbance  (MSD) 
transonic  theory.  The  wing  and  the  canard  are 
treated  in  a  sheared  fine  grid  system  that  is  embed- 
ed  in  a  global  Cartesian  crude  grid.  An  appropriate 
far  field  and  asymptotic  expression  for  the  velocity 
potential  derived  using  Green's  theorem  is  Imple¬ 
mented.  Results  are  presented  for  a  few  canard-wing 
configurations  and  compared  with  available  experi¬ 
mental  data.  The  weakening  of  the  wing  shock  due  to 
the  presence  of  the  canard  downwash  is  illustrated 
in  terms  of  contour  plots.  An  empirical  Incidence 
correction  for  the  wing  leading  edge  vortex  gives 
good  agreement  with  experiment  at  low  incidences. 

For  higher  angles  of  attack,  the  results  indicate 
that  a  more  sophisticated  vortex  roll-up  and  Induc¬ 
tion  model  is  required. 

Introduction 

A  number  of  highly  maneuverable  fighter  config¬ 
urations  such  as  the  H1MAT,  XFV-12A  and  Forward 
Swept  Hlng,  have  been  proposed  with  closely  coupled 
canard  systems  which  can  lead  to  several  advantages 
such  as  higher  trimmed-llft  capability.  Improved 
pitching  moment  characteristics,  and  reduced  trim 
drag.  Also,  the  geometric  characteristics  of  close¬ 
ly  coupled  canard  configurations  offer  an  improved 
longitudinal  cross-sectional  area  progression  which 
could  result  in  reduced  wave  drag  at  low  supersonic 
speeds.  Additional  advantages  associated  with  Im¬ 
proved  side  force  capability  have  been  reported  by 
Re  and  Capone[l].  The  associated  Interaction  due  to 
such  surfaces  with  the  wing  as  well  as  those  from 
conventional  tall  planes  involves  Important  non¬ 
linear  phenomena  in  the  supercritical  speed  regime. 
These  effects  can  significantly  change  spanwlse  load 
distributions  as  well  as  the  effective  Incidence 
field.  Corresponding  modifications  of  aerodynamic 
performance  and  stability  characteristics  are  there¬ 
fore  to  be  anticipated  not  only  for  fighter  config¬ 
urations,  but  with  tails  interacting  with  large- 
aspect-ratio  wings  typical  of  transport  arrange¬ 
ments  as  well. 


To  understand  clearly  the  closely  coupled 
canard-wing  Interference  transonic  flow  fields  with 
emphasis  on  benefits  associated  with  longitudinal 
and  lateral  positioning  that  have  been  discussed  In 
the  subsonic  regime  by  Lacey£2]  and  transonic  regime 
by  Gloss  and  Washbum{3],  It  Is  essential  to  develop 
a  computational  model  that  properly  simulates  the 
full  nonlinear  coupling  of  both  surfaces.  In  this 
respect,  the  mutual  Induction  arises  naturally  In 
the  computational  model  In  contrast  to  other 
"Inversion"  techniques  Mich  have  been  advocated. 

The  latter  employs  an  Iterative  procedure  In  which 
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each  surface's  isolated  aerodynamics  Is  alternately 
corrected  for  induction  from  the  other.  This  con¬ 
cept  may  have  some  validity  for  large  stagger  and 
gap  cases,  but  may  not  converge  for  closely  coupled 
arrangements.  In  fact,  even  If  convergence  1$  ob¬ 
tained,  the  result  may  not  be  the  physically  correct 
solution. 

A  mutually  interacting  two-dlmenstlonal  canard- 
wing  model  based  on  the  Karman-Guderley  transonic 
small  disturbance  theory  Is  reported  In  Ref.  [4] 
showing  several  parametric  studies  Involving  canard¬ 
wing  angles  of  attack,  stagger  and  oap  between  the 
two  lifting  elements.  The  objective  of  the  present 
study  Is  to  extend  the  two-dimensional  analysis  to 
three-dimensions  by  modifying  an  existing  small  dis¬ 
turbance  theory  transonic  code  for  wing-body  ccmbl- 
nations[5]  to  Include  the  presence  of  a  canard  or 
tail.  The  approach  Is  to  Introduce  a  separate 
sheared  fine  grid  box  for  the  canard  similar  to  the 
one  used  for  the  wing.  However,  the  experience  from 
the  two-dimensional  analys1s[4]  Indicates  that 
erroneous  solutions  can  result  from  the  use  of  two 
fine  grid  box  arrangement  for  very  closely  coupled 
canard-wing  systems  due  to  placement  of  the  outer 
boundaries  of  the  fine  grid  box  in  strong  nonlinear 
regions  of  the  other  surface.  To  avoid  this  situa¬ 
tion,  when  the  canard  Is  In  the  close  proximity  of 
the  wing,  a  single  fine  grid  box  Is  used  that  en¬ 
closes  both  the  lifting  surfaces.  By  contrast,  the 
two  fine  grid  box  arrangement  is  used  for  mildly 
coupled  canard-wing  configurations.  Beside  the  i>um- 
erlcal  Implementation,  the  paper  also  Includes  a 
detailed  analysis  of  the  appropriate  far  field  pro¬ 
cedure  based  on  an  asymptotic  theory.  The  present 
paper  neglects  the  downward  deflection  of  the  wake, 
vortex  roll-up  and  leading  edge  separation  phenom¬ 
ena.  These  can  be  Important  effects  and  are  being 
considered  In  separate  analyses. 

Results  are  shown  for  a  few  sample  canard-wing 
configurations  and  compared  with  available  experi¬ 
mental  data. 

Equations  am  Boundary  Conditions 

The  modified  form  of  the  transonic  classical 
small  disturbance  theory  In  conservation  form  as 
used  In  Ref.  [5],  will  be  eaploycd  In  the  current 
formulation  of  the  canard/tall  problem. 

[o  -  ^  nl7‘.;  ♦Mo;]) 

I  (1) 
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The  Move  equation  is  In  Me  conservation  form 
usually  preferred  for  bettor  shack  resolution.  Use 
of  Eg.  (1)  allows  one  to  apply  the  wing/ canard  sur¬ 
face  boundary  conditions  at  their  respective  assn 
chord  planes  and  transfer  boundary  cendltteni  on  the 
fuselage  to  a  convenient  computational  prl mastic 
surface.  Marring  to  Figure  1,  the  boundary  con¬ 
ditions  an  the  wing  and  canard  anna  plane  era  given 
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(2) 


where  z  -  zH|Ct  defines  the  wing/canard  mh  plane 

with  a  ♦  sign  denoting  the  upper  surface  and  a 
-  sign  the  lower  surface.  The  wing/canard  shape  is 
prescribed  by  c(x,y).  Application  of  boundary 

conditions  on  the  Mean  plane  avoids  the  need  for 
body- fitted  napping  procedures.  The  subscripts  w 
and  c  denote  the  wing  and  canard,  respectively. 

Across  the  vortex  wake,  there  will  be  a  Jump  In 
the  velocity  potential  d  which  will  be  Incorporated 
In  the  calculation  as  a  boundary  condition  while 
solving  Eq.  (1)  at  grid  points  neighboring  the  wake. 
The  wake  juwp  condition  Is 

[♦]HIN6  •  ♦(x,y,zw+)  -  d(xof,zw-)  \ 
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Fig.  1.  General  canard-wing-fuselage  arrangement 
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where  xTE  prescribes  the  trailing  edge  shape  and 
rw(y)  and  rc(y)  denote  the  sectional  circulation 
given  by 

rw(y)  -  ♦(*TE,wM,y,V) 

*  *(xTE,w<y)’y’V) 

* 

rc(y)  -  ♦(xTE<c(y).y.zc*) 

*  ♦<xTE.c<y),y,V) 


Fig.  2.  Coeputatlonal  domain  In  for  the  canard¬ 
wing-fuselage  arrangement 


In  general,  canard-wing  configurations  will  have 
different  planfona  shapes  for  both  lifting  units, 
and  In  particular,  different  leading  and  trailing 
edge  sweeps.  In  order  to  treat  arbitrary  arrange¬ 
ments,  It  Is  essential  to  Incorporate  separate 
shearing  transformations  for  both  surfaces  to  un¬ 
sweep  them  to  rectangular  planforms  In  a  computa¬ 
tional  domain.  This  can  be  accomplished  using  the 
following  shearing  transformation. 


^•e  "TE.w.c^  -  \E,w,c< 

Vc"y 
<w.c  "  1 


Since  the  transformation  far  the  wing  Is  different 
from  the  canard,  the  computational  he*  containing 
the  canard  will  be  different  from  that  containing 
the  wing.  This  Is  Illustrated  in  Figure  t. 

in  crier  to  establish  a  cenewntcatlen  tatwean 
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from  one  to  the  other,  a  global  crude  grid  will  be 


Fig.  3.  single  fine  grid  arr angaaent 


set  up  which  cm 
Figure  t,  fhogl 
Cartesian  system 
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Eq.  5.  To  achieve  a  desired  flow  resolution  for  the 
two  units,  a  very  fine  grid  with  clustering  near  the 
leading  and  trailing  edges  will  be  used  In  the 
canard  and  wing  boxes.  The  Cartesian  crude  grid 
will  extend  several  chord  lengths  away  fro*  the  con¬ 
figuration  where  far  field  boundary  conditions  can 
be  set.  For  very  closely  coupled  canard-wing 
systems  the  concept  of  two  fine  grid  box  arrangement 
can  lead  to  erroneous  solutions  as  well  as  numerical 
stability  problems  due  to  placement  of  the  outer 
boundaries  of  the  fine  grid  boxes  in  strong  non¬ 
linear  regions  of  the  other  surface.  This  problem 
was  experienced  in  the  earlier  two-dimensional  anal - 
ysis[4J.  Thus,  when  a  closely  coupled  canard-wing 
system  is  analyzed,  a  single  fine  grid  box  that  en¬ 
closes  both  the  canard  and  the  wing  Is  used  (the 
single  fine  grid  box  has  a  shearing  transformation 
based  on  the  wing  planform).  This  is  Illustrated 
in  Figure  3.  In  this  situation,  the  location  of 
the  canard  leading  and  trailing  edges  within  the 
wing  sheared  fine  grid  system  is  obtained  by  inter¬ 
polation.  Also,  the  fine  grid  may  have  to  be  speci¬ 
fied  manually  to  achieve  desired  grid  clustering  for 
the  canard.  The  double  fine  grid  box  arrangement  of 
Figure  2  should  work  very  well  for  mildly  coupled 
canard-wing  systems. 

Solution  Procedure 

Equation  (1)  in  its  Cartesian  form  is  solved  in 
the  global  crude  grid  while  the  transformed  equation 
in  terms  of  the  new  variable  c.n.C  is  solved  in  the 
fine  grid  box  using  a  successive  line  over  relaxa¬ 
tion  algorithm.  After  each  relaxation  cycle  along 
the  crude  grid-fine  grid  interface,  the  potential 
values  from  the  crude  grid  are  Interpolated  on  to 
the  fine  grid  points  and  used  as  Olrlchlet  boundary 
condition  for  the  fine  grid  box.  Similarly,  at  the 
wing  and  canard  mean  planes,  the  potential  values 
for  the  crude  grid  points  are  obtained  from  the  grid 
values  and  provide  a  Olrlchlet  boundary  condition 
for  the  crude  grid  domain.  The  use  of  crude  grid- 
fine  grid  successive  sweep  procedure  accelerates  the 
convergence  process(6]. 

Far  Field  Solution 

Along  the  outer  boundaries  of  the  computational 
domain  It  is  Important  to  use  appropriate  far  field 
boundary  conditions.  These  will  be  obtained  from  a 
bi-wing  generalization  of  the  asymptotic  analysis  of 
Human  and  Cole[7]  for  non-lifting  airfoils  and  its 
three-dimensional  wing  extension  by  Klunker[8]. 


Fig.  4a.  Canard-wing  system  and  far-field 
boundaries 


Referring  to  Figure  4a,  a  bi-wing  system  without 
a  fuselage  is  shown.  The  far  field  derivation  pro¬ 
cedure  will  be  indicated  for  these  elements,  and  an 
extension  to  inclusion  of  the  fuselage  will  be 
briefly  Indicated.  Shown  in  the  ftgure  Is  the 
canard  surface,  Sc,  upstream  of  the  main  wing,  S^. 

Also  depicted  are  the  respective  trailing  vortex 
sheets,  $y^  and  $y,  of  bot*“  elements.  Ai*0  pro¬ 
jection  of  the  compound  lifting  system  is  indicated 
in  Figure  4b. 

For  this  derivation,  we  neglect  the  dowm.,rd 
deflection  of  the  wake  for  moderately  large  down¬ 
stream  locations  of  the  computational  approximation 
uf  the  Trefftz  plane  boundary.  Also  neglected  in 
this  treatment  are  vortex  roll-up  and  leading  edge 
separation  phenomena. 

Using  the  procedure  from  Refs.  [7]  and  [8],  we 
obtain  the  aforementioned  integrodlfferential  equa¬ 
tion  for  the  velocity  potential  4  by  applying 
Green's  theorem  to  the  small  disturbance  equation. 
Here  we  deal  with  the  Karman-Guderley  formulation 
to  illustrate  the  bas^c  method  to  be  applied  to  a 
modified  small  disturbance  solver[4].  Using  suit¬ 
able  strained  coordinates  (x,y,z)  Indicated  in  the 
figure  (for  field  points  referred  to  an  origin  at 
the  center  of  the  wing  leading  edge,  and  a  corre¬ 
sponding  set  U.n.c)  for  the  dunay  variables,  the 
small  disturbance  equation  for  the  perturbation 
potential  4  in  this  system  is 


For  purposes  of  application  of  Green's  theorem,  we 
introduce  a  fundamental  solution  G  which  satisfies 
the  following  equation 

Vt.n.cG  *  6K-x)«(n-y)6(c-z)  •  (7) 

Thus,  G  is  the  usual  unit  source  given  by 

6  "  '  A 

R*  “  (x-5)1  +  (y-n)a  ♦  (z-c )*  • 

Using  Green's  theorem  on  the  domain  cut  along 
the  approximate  wake  locations  shown  in  Figure  4b, 
i.e.,  the  planes,  y ■  0  and  y  ■  0,  we  obtain 
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where,  In  the  partial  integrations,  all  terms  from 
the  shock  surfaces  such  as  o  (shown  schematically  in 
Figure  4a)  vanish,  and 

[f(K,y.z)j  =  f(x.y, 0*)  -  f(x,y,0-) 

and  the  C.c  subscripts  denote  partial  differentia¬ 
tion  In  the  Integrals. 

The  Integrals  over  the  finite  doaalns  can  be 
simplified  by  approximating  the  kernels  for  x2  +  y2 
♦  z1-*-,  with  normalized  semispan  B  ■  o(y2  +  z2). 
This  approximation  and  the  condition  that  can  be 
arbitrarily  prescribed  leads  to  the  desired  asymp¬ 
totic  formula 

♦  ■  IL  ♦  IT  ♦  lfC>  ♦  l{c)  ♦  Iv  (9) 

where.  If  the  wing  and  canard  are  given  by  the 
equations, 

1  ’  fy  j(*«y)  ,  ^  ^  ,  -B  <  y  <  B) 

2  “  f^i(x,y)  •  (xle><x<xTE)  *  •Bc<y<Bc) 
respectively,  then 
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and  B,  Bc  representing  respectively  the  scaled  wing 
and  canard  semi spans. 

Note  that  In  (9) 
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and  the  rest  of  the  terms  are  higher  order  being 
0(Rq  ),  which  can  be  neglected  for  sufficiently 

large  computational  domains.  It  is  therefore  appar¬ 
ent  that  owing  to  the  linearity  of  the  far  field, 
the  dominant  effect  Is  the  sum  of  the  two  lift  con¬ 
tributions  IL  and  Ij.  Where,  as  In  the  two- 

dimensional  case  discussed  in  Ref.  [4],  there  Is 
however  nonlinear  coupling  between  the  circulations 
r  and  rc. 

In  a  similar  fashion,  the  effect  of  the  for 
field  of  a  slender  body  contribution  can  be  added  to 
(9)  using  the  nonlinear  line  source  far  field 
expression  given  In  Ref.  [9].  Denoting  this  term  as 
Ip,  we  have  that 

Ip  «  o(Sq/Rq)  as  Rq  ♦  ■ 

where  SQ  is  the  body  base  area.  If,  as  In  typical 
cases,  Sq/Sc  «  1  and  Sc  «  S,  then  Ip  will  be 
neglected. 

In  the  vicinity  of  the  projections  of  the  lift¬ 
ing  elements  in  the  Trefftz  plane,  the  previous 
approximations  become  Invalid.  The  appropriate 
'limits  are  obtained  by  assuming 

z»y »zc  -  0(B,Bc)  as  x  ♦  •  (11) 

Using  the  limit  (11),  (8)  can  be  approximated  to 
give  the  contribution  of  the  projection  of  the  vor¬ 
tex  sheets  from  both  elements 
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for  the  fuselage  absent. 

Equation  (12)  corresponds  to  the  solution  of  the 
cross  flow  Laplace's  equation  for  the  perturbation 
potential  In  the  Trefftz  plane  cut  about  the  projec¬ 
tions  of  the  Ming  and  canard.  If  a  fuselage  Is 
present,  the  solution  can  be  obtained  by  superposing 
the  appropriate  Trefftz  plane  Unit  of  the  nonlinear 
line  source  described  In  Ref.  [9]. 

Results 


the  canard.  To  natch  the  experimental  lift,  the 
conputer  code  was  run  at  an  angle  of  attack  of  5.8*. 
while  the  experlnent  was  conducted  at  4.18*.  The 
Increased  angle  of  attack  used  In  the  conputatlonal 
run,  Is  an  ad  hoc  fix  to  account  for  the  suction 
peak  associated  with  the  weak  leading  edge  vortex 
present  In  the  experlnental  run  as  well  as  viscous 
and  wall  interference  effects.  Still,  there  Is  a 
discrepancy  In  the  surface  pressure  distribution 
conparlson  between  nunerlcal  prediction  and  experl¬ 
nental  data^as  shown  In  Figure  7.  The  presence  of 
the  canard  produces  a  downwash  field  in  front  of  the 
wing  causing ^an  effective  decrease  In  the  angle  of 
attack.  This/results  in  a  reduction  In  the  nose 
pressure  distribution.  The  conparlson  on  the  lower 
surface  pressure  with  and  without  the  canard  Is  very 
9ood.  At  higher  angles  of  attack  the  conparlson 
with  data  Is  very  poor  due  to  the  Inability  of  the 
conputatlonal  nodel  to  predict  the  leading  vortex 
systen.  The  case  shown  In  Figure  7  at  least  Indi¬ 
cates  the  right  trend  predicted  by  the  conputatlonal 
model  with  the  canard  present  and  Interestingly,  Is 
In  qualitative  agreement  with  the  ltehavlor  of  two- 
dimensional  systems  treated  In  Ref.  [4]. 


All  the  canard-wing  calculations  presented  here 
were  performed  using  the  LBL  CDC  7600  machine.  A 
single  fine  grid  encompassing  both  the  canard  and 
wing  was  used  since  the  configurations  consider¬ 
ed  in  this  paper  were  all  very  closely  coupled 
canard  systems.  A  typical  fine  grid  consisted  of 
90x30x20  points  and  the  crude  grid  30x20x20  points. 

A  total  of  100  crude-fine  cycles  requiring  10-15 
minutes  of  CPU  time  usually  resulted  in  satisfactory 
convergence  in  which  the  difference  In  $  between  two 
successive  cycles  was  of  the  order  of  10  . 

Gloss  and  Washburn[3]  have  made  extensive  exper¬ 
imental  measurements  on  the  close  coupled  canard- 
wing  configuration  shown  In  Figure  5.  The  wing 
and  canard  surfaces  consisted  of  sharp  leading  edge 
circular  arc  airfoils.  The  pressure  data  In  their 
report  clearly  indicate  a  strong  leading  edge  vortex 
system  especially  at  higher  angles  of  attack.  Since 
the  computational  model  In  its  present  form  cannot 
predict  leading  edge  vortex  flows,  for  conparlson 
purposes,  only  the  small  angle  of  attack  case  Is 
considered,  where  the  leading  edge  vortex  is  weak. 
Figure  6  shows  the  comparison  of  spanload  distribu¬ 
tion  on  the  wing  with  and  without  the  presence  of 


Another  canard-wing  configuration  considered  Is 
shown  In  Figure  8.  Two  different  canard  positions 
In  the  vertical  direction  are  analyzed.  For  this 
case,  the  effect  of  the  canard  and  Its  position  on 
the  wing  sectional  lift  characteristics  Is  shown  In 
Figure  9,  for  «  0.84  and  angle  of  attack  of  3*. 
The  lift  produced  by  the  canard  produces  a  downwash 
in  front  of  the  wing  causing  the  wing  lift  to 
decrease.  However,  at  wing  span  stations  away  from 
the  canard  tip,  the  effect  of  the  canard  Is  minimal . 
Corresponding  to  this  case,  the  pressure  contours  at 
the  wing  span  station  (y/by)  -  0.245  Is  shown  in 
Figure  10  for  both  the  canard  off  and  on  situations. 
Figure  10(a)  shows  the  Isolated  wing  case  (canard 
off)  exhibiting  a  shock  on  the  upper  surface.  When 
the  canard  Is  placed  at  0.5  Cp  above  the  wing  plane 
(Cp  Is  the  wing  root  chord),  the  shock  weakening  on 
the  wing  due  to  downwash  produced  by  the  canard  Is 
clear  In  Figure  10(b).  Even  further  weakening  of 
the  shock  due  to  closer  positioning  of  the  canard  Is 
seen  In  Figure  10(c).  Figure  11  shows  the  surface 
pressure  distribution  on  the  wing  at  various  span 
stations  with  canard  off  and  on.  This  figure  also 
clearly  shows  the  reduction  In  wing  shock  strength 
when  the  canard  Is  present. 
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Conclusion! 

A  simple  canard-wing  transonic  cosputational 
modal  Is  presented.  The  formulation  neglects  wake 
deflection,  vortex  roll-up  and  leading  edge  separa¬ 
tion  phenomena.  An  appropriate  far  field  expres¬ 
sion  for  the  velocity  potential  is  derived  and  im¬ 
plemented.  Results  are  presented  for  txo  canard- 
wing  configurations  and  some  conparison  with  exper¬ 
imental  data  is  made.  An  meplrical  Incidence 
correction  for  the  wing  leading  edge  vortex  gives 
good  agreement  with  experiment  it  low  incidences. 

The  presence  of  the  canard  produces  a  downwesh  field 
that  tends  to  weaken  the  shock  system  on  the  wing 
surface.  For  higher  Incidence,  soma  kind  of  a  lead¬ 
ing  edge  potential  vortex  core  nodal  superimposed 
onto  the  present  computational  treatment  would 
greatly  enhance  the  prediction  capability.  This 
will  be  pursued  in  the  future  along  with  suitable 
models  for  Incorporating  the  wake  roll-up  phenomena. 
The  computational  speed  can  also  be  enhanced  by  use 
of  Implicit  factored  algorithms. 
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Fig.  10.  Constant  pressure  contours  at  span  station 
(y/bw)  «  0.245.  =  0.84,  c^>e  -  3° 


Fig.  11.  Surface  pressure  distribution  on  the  wing 
with  canard  off  and  on.  IL  ■  0.84, 


